Abstract. A detailed study of the planetary nebula 
NGC 6565 has been carried out on long-slit echellograms 
(A/AA=60000, spectral rangc=AA3900-7750A) at six, 
equally spaced position angles. The expansion velocity 
field, the c(H/3) distribution and the radial profile of the 
physical conditions (electron temperature and density) are 
obtained. The distance, radius, mass and filling factor of 
the nebula and the temperature and luminosity of the cen- 
tral star are derived. The radial ionization structure is 
analyzed using both the classical method and the photo- 
. ionization code CLOUDY. Moreover, we present the spa- 
^? ' tial structure in a series of images from different directions, 
^1;;^ I allowing the reader to "see" the nebula in 3-D. NGC 6565 
■ results to be a young (2000-2500 years), patchy, optically 
thick triaxial ellipsoid (a=10.1 arcsec, a/b=1.4, a/c=1.7) 
projected almost pole-on. The matter close to major axis 
: was swept-up by some accelerating agent (fast wind? ion- 
ization? magnetic fields?), forming two faint and asym- 
, metric polar cups. A large cocoon of almost neutral gas 
completely embeds the ionized nebula. NGC 6565 is in a 
recombination phase, because of the luminosity drop of 
the massive powering star, which is reaching the white 
dwarf domain (logT* -5.08 K; logL^/Lg -2.0). The stel- 
lar decline started about 1000 years ago, but the main 
: nebula remained optically thin for other 600 years before 
' the recombination phase occurred. In the near future the 
ionization front will re-grow, since the dilution factor due 



> 



o 

' to the expansion will prevail on the slower and slower stel- 



lar decline. NGC 6565 is at a distance of 2.0(±0.5) Kpc 



O and can be divided into three radial zones: the "fully ion- 



ized" one, extending up to 0.029-0.035 pc at the equator 



^ ' (0.050 pc at the poles), the "transition" one, up to 0.048- 
'. 0.054 pc (0.080 pc), the "halo", detectable up to 0.110 
' P*^- The ionized mass (~0.03 Mq) is only a fraction of the 
, total mass (> 0.15 Mq), which has been ejected by an 
equatorial enhanced superwind of 4(±2)xl0^'^ Mq yr~^ 
. 5^ ^ lasted for 4(±2)xl03 years. 
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1. Introduction 

The late evolution of low and intermediate mass stars (1.0 
Mq <M* <8.0 Mq) is characterized by the planetary 
nebula (PN) metamorphosis: an asymptotic giant branch 
(AGB) star gently pushes out the surface layers in a flori- 
legium of forms, crosses the HR diagram and reaches the 
white dwarf regime (Aller 1984, Pottasch 1984, Osterbrock 
1989). 

So far the interpretation of the exuberant morpholo- 
gies exhibited by the PNe in terms of detailed three- 
dimensional structures and physical conditions of the ion- 
ized gas was limited by projection effects, leading to ap- 
proximate spatial forms and to unrealistic assumptions for 
the main parameters, e.g. electron temperature and elec- 
tron density constant all over the nebula (Aller 1984, 1990, 
1994). 

A 3-D reconstruction technique for studying at large 
and small scales the morphology, physical conditions, ion- 
ization, spatial structure and evolutionary status of PNe 
has been introduced by Sabbadin et al. (2000a, b) and 
Ragazzoni et al. (2001), based on echellograms of moder- 
ate spectral resolution (R~ 22000 - 25000). The key of 
the 3-D methodology is simple: the PN is an expanding 
plasma. Thus the position, thickness and density of each 
elementary volume can be derived from the radial velocity, 
width and flux of the corresponding emission. 

The procedure has been applied to: 

- NGC 40: an optically thick, very low excitation barrel- 
shaped nebula with thin arcs emerging at both ends 
of the major axis, powered by a luminous and "cold" 
^^T<S star i)i (^s<-nting a large mass-loss rate. The fast. 
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Fig. 2. [NII]/[OIII] distribution over NGC 6565 (original 
frames and orientation as in Fig. 1), showing the large 
stratification effects and the faint, low ionization regions 
protruding from the main nebula in PA~145°. 

hydrogen depleted photosphcric material ejected by 
the nucleus is gradually modifying the chemical com- 
position of the innermost nebular regions (Sabbadin et 
al. 2000a); 

- NGC 1501: an evolved, high excitation, optically thin 
oblate ellipsoid, denser in the equatorial belt, deformed 
by several bumps, embedded in a homogeneous, in- 
wards extended cocoon and ionized by a " hot" and 
luminous WC4 star exhibiting nonradial g-mode pul- 
sations (Sabbadin et al. 2000b; Ragazzoni et al. 2001). 

In order to deepen the analysis, we have started a 
survey at high spectral and spatial resolutions with the 
ESO NTT. The superb quality of this material allows us 
to study at unprecedented accuracies objects with differ- 
ent morphology, e.g. NGC 7009 (the "Saturn" nebula), 
the tetra-lobed IC 4634, the butterfly HB 5, Mz 3 and 
NGC 6537, the double-envelope NGC 5882, NGC 6153 
and NGC 6818. 

To simplify the application of the 3-D method we de- 
cided to begin with an "easy" nebula without FLIERS 
(fast, low ionization emitting regions), BRETS (bipolar, 
rotating, episodic jets), ansae, wings, multiple envelopes 
etc., and naively selected NGC 6565. 

2. The nebula 

The bright, compact PN NGC 6565 (PNG 003.5-04.6, 
Acker et al. 1992) is "a minute oval ring 10"x8" in p. a. 
about 5°. Considerably fainter along the major axis, and 
the center is relatively vacant" (Curtis 1918). Its HST ap- 
pearance is shown in Fig. 1: the [OIII] emission forms a 
quite homogeneous oval ring, sprinkled by a number of 
dark globules and spikes, mainly in the Northern sector. 
It is surrounded by a sharp, patchy [Nil] layer present- 
ing two faint lobes in PA~145°. Note also the light N-S 
asymmetry of the nebula with respect to the (barely vis- 
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Fig. 1. WFPC2 appearance of NGC 6565 in [OIII] (left) and [Nil] (right). The images have been retrieved from the 
HST archive. North is up and East is to the left. 



ible) central star. The large stratification effects and the 
external, faint, low ionization regions are better seen in 
Fig. 2, presenting the [Nil] /[OIII] distribution. 
NGC 6565 appears peculiar in many respects: 

- despite the simple elliptical shape (Stanghellini et al. 

1993; Gorny et al. 1997), it shows a rich low ionization 
spectrum typical of bipolar PNc (Greig 1972) and is N 
overabundant (Aller et al. 1988); 

- the infrared and radio excesses suggest the presence of a 

large amount of dusty neutral gas causing a detectable 
absorption of the optical radiation (Gathier et al. 1986; 
Stasinska et al. 1992); 

- in spite of the high surface brightness of the nebula, 

indicative of an early evolutionary phase, the central 
star is faint and hot, i.e. it is in, or close to, the white 
dwarf domain. Following Tylenda (1986), NGC 6565 
could be in a recombination phase, due to the recent 
luminosity drop of the powering star. 

Although the nebula is sometimes listed among the 
distance calibrators (Pottasch 1984; Sabbadin 1986; Gath- 
ier 1987; Zhang 1995), its distance is still poorly defined 
and the distance-extinction law gives contradictory results 
(Gathier et al. 1986, Maciel et al. 1986 and de Oliveira- 
Abans & Faundez-Abans 1991). 

This paper is structured as follows: Sect. 3 presents the 
observational material and the reduction procedure. Sect. 
4 is dedicated to the gas kinematics. Sect. 5 contains the 
Ha/H/3 spectral maps, Sect. 6 analyses the radial distribu- 
tion of the physical conditions (electron temperature and 
electron density), in Sect. 7 we derive the distance, the 
nebular mass and the central star parameters, in Sect. 8 
the radial ionization structure and the distribution of the 
gas (ionized and neutral) are presented. Sect. 9 concerns 
the application of the photo- ionization model (CLOUDY), 
Sect. 10 describes the 3-D structure of the nebula in dif- 
ferent ions, a short discussion and the conclusions are pre- 
sented in Sect. 11. 

3. The observational material and reduction 
technique 

Six long-slit echellograms of NGC 6565 (exposure time 

900s; spectral resolution 60000 with a slit 1.0 arcsec wide) 
have been obtained on July 28, 2000 at the ESO NTT 
equipped with EMMI in REMD mode using grating #14, 
grism #3 and CCD #36 (spectral range AA3900-7750A, 
mean spatial scale along the slit 0.27 arcsec pix~^). During 
the observations the sky transparency was excellent and 
the seeing fluctuated between 0.53 and 0.89 arcsec. The 
30 arcsec long slit was centered on the nebula for all the 



selected PA, ranging from 0° to 150° with a constant step 
of 30°. 

The reduction method follows conceptually the stan- 
dard procedure, including bias, flat field, distortion cor- 
rection and wavelength and flux calibration. 

The slit length, much longer than the order separation, 
leads to a considerable superposition of the echelle orders. 
This, on the one hand does not affect the emission line 
spectrum of NGC 6565, on the other hand it does the 
flat fleld exposures. Thus, we were forced to apply a zero- 
order correction using a white-light flat field. Thanks to 
the excellent CCD uniformity, the pixel to pixel sensitivity 
fluctuation left by this procedure is below 5% (verifled 
on the night sky and the comparison spectrum lines). All 
things considered, it is a modest price to pay for keeping 
the whole spatial information of the PN emissions in the 
AA3900-7750A spectral range. 

After bias and flat field correction, the orders con- 
taining one or more emissions of physical interest were 
extracted and treated as single, long- slit spectra. The 
distortion correction and the wavelength calibration were 
performed by means of Th-Ar lamp exposures. During the 
frame registration, the various orders were scaled to the 
slit height at Ha, chosen as reference line. Echellograms 
of the spectrophotometric standard star HR5501 (sampled 
at several positions along each order) were used to obtain 
the flux calibration. 

We stress the substantial difference between our obser- 
vational procedure and the one generally adopted, which 
introduces an interference filter to isolate a single order. 
On the contrary, we cover the whole AA3900-7750A spec- 
tral range with a single exposure. 

For illustrative purposes, in Fig. 3 we present the main 
part of the NTT-hEMMI frame of NGC 6565 at PA=150°, 
which includes the external, low ionization regions, and in 
Fig. 4 the detailed structure (at the same PA) of A4686A 
of Hell, A5007A of [OIII], A5876A of Hel, A6300A of 
[01], A6563A of HI, A6584A of [Nil], A673lA of [SII] and 
A7135A of [Arlll]. They suggest that the faint, low ion- 
ization regions visible in Fig. 2 are the projection of ex- 
tended cups protruding from the main nebula. Note also 
the large stratification effects, the [01] rises in the out- 
ermost parts and, in spite of the blurred Ha appearance 
due to thermal motions plus expansion velocity gradient 
plus fine structure, the similarity of the HI and [OIII] dis- 
tributions. Following Williams (1973), these features are 
indicative of an optically thick nebula powered by a high 
temperature star. 

A first, qualitative picture of the NGC 6565 struc- 
ture can already be inferred from Figs. 1 to 4: the ob- 
ject is an ellipsoid with extended polar cups, projected 
almost pole-on. The possible optical appearance of the 
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Fig. 3. Main part of tlie NTT+EMMI frame of NGC 6565 at PA=150° (in logarithmic scale); the principal nebular 
emissions are identified and the strongest night sky lines (N.S.) indicated. The spatial scale is given by the N.S. 
extension, corresponding to the slit length of 30". 

Fig. 4. Detailed structure (in logarithmic scale) of some representative emissions observed in NGC 6565 at PA=150°. 
The orientation is as in Fig. 3. The fluxes are multiplied for the factor given in parenthesis to render each emission 
comparable with A5007A of [OIII]. The A6300.304A night sky line was removed from the [01] nebular emission. Part 
of the [01] A6363A emission (cchoUc order 96) appears as a moustache in the lower-left corner of the A6300A image; 
(echelle order 97). The fast, blue-shifted wisp in the lower-left corner of the [Nil] and [SII] frames is real. A hot pixel 
column grazes the Ha emission. The zero- velocity pixel column (zvpc) and the central star pixel line (cspl), defined 
in Sect. 4, are shown in the lower-right panel. 



nebula seen equatorial-on is illustrated in Fig. 5, show- 
ing the PN NGC 6886 in [OIII] and [Nil]. NGC 6565 and 
NGC 6886 are almost twins in many respects: nebular 
spectrum, physical conditions, ionization structure, evolu- 
tionary phase, luminosity and temperature of the central 
star. Likely, they derive from similar progenitors. 

4. The gas kinematics 

We introduce two important definitions relative to each 
emission line, which will be extensively used in the paper 
(they are shown in the lower-right panel of Fig. 4): 

1 ) the zcro-vclocity pixel column (zvpc, see Sabbadin ct 

al. 2000a, 2000b, and Ragazzoni et al. 2001), corre- 
sponding to the recession velocity of the whole nebula, 
collects the photons of the ionized gas expanding per- 
pendicularly to the line of sight; 

2 ) the central star pixel line (cspl), parallel to the disper- 

sion, represents the nebular material projected at the 
apparent position of the star, whose motion is purely 
radial. 

In a certain sense the zvpc and the cspl are comple- 
mentary: the first is independent on the expansion velocity 
and gives only spatial information, the second measures 
the velocity field and contains only kinematical informa- 
tion. By combining the zvpc and the cspl we will derive 
the overall spatio-kinematical properties of the nebula. 

The expansion velocity (Vexp) of the ionized gas in 
NGC 6565 was derived from the analysis of the cspl in 
the different ions. Note that the cspl is the same in all 
the frames, the spectra being radially arranged. Thus, in 
order to achieve the faintest emissions (like A5198-5200A 
of [NI], A5517-5537A of [ClIII] and A7005A of [ArV]), the 
six echellograms were combined. 

The results are contained in Table 1, where the ions arc 
put in order of increasing ionization potential (LP-). Be- 
sides the peak separation (2Vexp) of the main nebula (col- 
umn 4), in column 5 we give Vexp of the faint, fast, blue- 
shifted component emitted by the approaching, southern 
polar cup, clearly visible in Figs. 3 and 4. 

As the cspl refers to the kinematical properties of the 
(almost) polar gas in NGC 6565, so the zvpc contains the 



spatial properties of the equatorial matter. The intensity 
peak separations, 2rzvpc, in the different emissions are pre- 
sented in columns 6 to 1 1 of Table 1 . 

Typical errors for 2Vexp(main) and 2rz;vpc are 1.0 km 
s~^ and 0.15 arcsec for the strongest forbidden emissions 
([01], [SII], [Nil] and [OIII]) to 3.0 km s^^ and 0.3 arcsec 
for the faintest ones ([NI], [Oil], [ClIII] and [ArV]). The 
corresponding inaccuracies for the recombination lines are: 
2.0 km s^^ and 0.2 arcsec for HI, 1.5 km s~^ and 0.2 arcsec 
for Hel and 3.0 km s^^ and 0.3 arcsec for Hell. The error 
in Vexp (cup) (column 5) is 3.0 km s~-^ in all cases. 

The correlation bcrtwccn cspl and zvpc, evident in Ta- 
ble 1, is summarized in Fig. 6, showing both Vexp and 
fzvpc vs. I. P.. For reasons of homogeneity, the peak sepa- 
rations at different PA are normalized to the [01] value. 

NGC 6565 follows the Wilson's law ("the high- 
excitation particles show smaller separations than the low 
excitation particles" , Wilson 1950) with some caveats: 

- the "high ionization" zone (I.P.>50 eV) is poorly de- 

fined, containing only two emissions: A7005A of [ArV] 
is exceedingly faint and A4686A of Hell consists of thir- 
teen fine structure components whose spectral distri- 
bution presents a noticeable blue-shifted tail lowering 
the measurement reliability (in fact, 2 Vexp (Hell) from 
A6560A is 38±3 km s"^); 

- the "medium ionization" region (20 eV<I.P.<50 eV) is 

characterized by a steep velocity gradient. Vexp and 
Tzvpc in [ArlV] are considerably lower than in the other 
ions, as fomid by Wilson (1950) for NGC 3242 and 
NGC 7662. On the contrary, [Nelll], a "medium- high" 
ionization species, is emitted in the "medium-low" ex- 
citation region, as observed in NGC 6818, NGC 6886 
and NGC 7662 by Wilson (1950), and in NGC 6720 
by Reay & Worswick (1977). These peculiarities will 
be analyzed in Sect. 9 (photo-ionization model); 

- the "low ionization" zone (I.P.<20 eV) is steep and well 

defined, with the obvious exception of hydrogen, which 
is emitted by the whole ionized nebula, independently 
on the plasma conditions. 

Assuming a minimum degree of symmetry within the 
nebula we can assess that the expansion velocity is pro- 
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Fig. 5. HST imaging of the PN NGC6886 in [OIII] and [Nil], showing the possible appearance of NGC 6565 seen 
almost equatorial-on. 



Table 1. Expansion velocity in the cspl and peak separation in the zvpc for NGC 6565. 



Ion 


A 

(A) 


LP. 

(eV) 


2Vcxp(niain) 


Vexp(cup) 






peak separation (") 






(km/s) 


(km/s) 


PA=0° 


PA=30° 


PA=60° 


PA=90° 


PA=120° 


PA=lf 


[NI] 


5198 


0.0 


70.9: 


- 


8.7: 


8.4: 


8.0: 


7.6: 


8.0: 


8.3: 


[01] 


6300 


0.0 


67.7 


55 


8.7 


8.4 


8.0 


7.6 


8.0 


8.3 


[SII] 


6731 


10.4 


65.0 


54 


8.5 


7.8 


7.5 


7.3 


7.5 


7.9 


[Oil] 


7319 


13.6 


59.4: 




7.8: 


7.7: 


7.2 


7.0 


7.4 


7.8: 


HI 


6563 


13.6 


55.8 


48: 


7.1 


6.7 


6.2 


6.0 


6.2 


6.5 


[Nil] 


6584 


14.5 


60.6 


54 


8.2 


7.6 


7.4 


7.1 


7.4 


7.8 


[SIII] 


6312 


23.4 


57.5 




7.6 


6.9 


6.7 


6.3 


6.5 


6.7 


[ClIII] 


5537 


23.8 


57.9: 




7.0: 


6.9: 


6.5: 


6.1: 


6.3: 


6.7: 


Hel 


5876 


24.6 


58.7 


52: 


7.5 


6.8 


6.5 


6.2 


6.5 


6.8 


[Arlll] 


7135 


27.6 


57.2 


50 


7.3 


7.0 


6.8 


6.3 


6.5 


7.0 


[OIII] 


5007 


35.1 


56.7 


49 


6.9 


6.4 


6.3 


6.0 


6.3 


6.4 


[ArlV] 


4740 


40.7 


32.4 




3.9 


3.6 


3.5 


3.5 


3.6 


3.9 


[Nelll] 


3967 


41.0 


55.2 




7.2 


6.6 


6.5 


6.3 


6.6 


6.6 


Hell 


4686 


54.4 


44.2: 


42: 


4.5 


3.5: 


3.8: 


3.7 


4.0 


4.1 


[ArV] 


7005 


59.8 


17.0: 






3.4: 


3.2: 


2.7: 


3.0: 


2.7: 



portional to the distance from the central star, in quali- 
tative agreement with the results of the classical papers 
of Wilson (1950) and Weedman (1968), performed on a 
consistent sample of objects. 

The exact value of the proportionality constant, 
A=Vexp/r", remains unknown, since at the moment no 
point of the nebula possesses both the expansion velocity 
and the distance from the star. Although a rigorous treat- 
ment of the problem is presented in Sect. |^, an intuitive 
solution is evident in Fig. 4: by chance the combination of 
spatial scale and spectral dispersion in the NTT-I-EMMI 
echellograms gives emission line structures of NGC 6565 
reproducing the tomographic maps in the slice covered by 
the slit. 

Let's come back to the gas kinematics. Further infor- 
mation is in the position-velocity maps, i.e. in the com- 
plete radial velocity field. The maps, shown in Fig. 7 for 
the six PA, are in the nebula rest frame. The systemic 
heliocentric velocity is Vr0 — —5.2 ± 0.5 km s~^. We 
have selected Hell, [OIII] and [Nil] as representative of 
high, medium and low ionization regions, respectively, and 
Ha, which gives the overall distribution of the ionized gas 
(F1ux(Hq;) oc Ne2ocN(H+)2, where Ne and N(H+) are the 
electron and proton densities, respectively). 

A number of interesting features can be identified from 
these maps: the complex expansion velocity field charac- 
terized by the polar cups, the fast low excitation blue- 
shifted wisp in PA=0° and 150°, the large stratification 
effects, the blurred appearance of Ha due to a combina- 
tion of thermal motions, fine-structure and expansion gra- 
dient, and the blue-shifted tail in Hell (asymmetric distri- 
bution of the thirteen fine-structure components). In ad- 
dition, our spectra point out the presence of an external. 



very faint (untilted?) emission having different kinemati- 
cal properties with respect to the main nebula. This "halo" 
(barely visible in Figs. 2, 3 and 7) is more evident in [Nil] 
than in Ha and [OIII]. It appears double-peaked in [Nil] 
(2Vexp=25±3 km s^^), with the blue-shifted component 
stronger than the red-shifted one in the southern half of 
the nebula, and vice-versa in the northern half. In both 
Ha and [OIII], a single, broad emission (FWHM=25±5 
km s^^) is visible. 

The existence of a faint, low ionization envelope around 
NGC 6565 supports the Tylenda's (1983, 1986) suggestion 
that the nebula is in a recombination phase. 

The only kinematical study of NGC 6565 reported in 
the literature is by Meatheringham et al. (1988), who de- 
rived 2Vexp=41.8 km s'^ in [OH], 28.8 km s'^ in [OIII] 
and 16.6 km s~^ in Hell from the peak separation, and 
2Vexp [OIII] =54.4 km s"^ from the width at 10% maxi- 
mum intensity (following Dopita et al. 1985, this gives to 
a good approximation the highest expansion velocity). 

The disagreement with the results of Table 1 can be 
ascribed to the different spatial and spectral accuracies 
(Meatheringham et al. 1988 used a slit width of 3.9" and a 
spectral resolution ~26000), and to the different reduction 
techniques. To be noticed that the width at 10% maximum 
intensity in our echellograms gives 2Vexp[OIII]=90(±4) 
km s~^. Further comments will be contained in subsection 
10.2. 

5. The Ha/H/3 spectral maps 

In general (see Aller 1984, Pottasch 1984, and Osterbrock 
1989), the PN extinction is obtained by comparing the 
observed Balmer decrement (in particular Ha/H/?) to the 
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Fig. 6. Polar expansion velocity (left) and equatorial peak separation (right) vs. LP. in NGC 6565, showing the linear 
relation between Vexp and distance from the central star. The intervals of the peak separations at different PA, 
normalized to the [01] values, are reported. The ions deviating from the normal (cfr. text) behavior are marked. 



Fig. 7. Position-velocity maps for representative ions of NGC 6565 at the six PA observed. The orientation of each 
PA is indicated in the Hell panel. The contour step is AlogF=0.30. 



intrinsic value given by Brocklehurst (1971) and Hummer 
& Storey (1987). The estimates of (HQ!/H/3)obs reported 
in the literature for NGC 6565 span the range 3.79 (de 
Freitas Pacheco et al. 1992) to 4.65 (Aller et al. 1988). 
They all represent mean values over the slice of nebula 
covered by the slit. 

The spatial and spectral accuracies achieved by the 
superposition technique used (±0.15 arcsec and ±1.0 km 
s~^, respectively) allow us to extend the Ha/H/3 analysis 
to the whole spectral image. Note that the fine structure 
of the lines (the spectral distribution and the relative in- 
tensities of the seven components of Ha differ from the 
corresponding quantities of the as many components of 
H/3, Clegg et al. 1999) is uninfluent for an expanding gas 
at Te~10^ K. 

Fig. 8 shows the Hq:/H/3 isophotal contours su- 
perimposed to the Ha spectral emission (recall that 
F1ux(Hq;)ocN(H"'")^) for some representative PA of 
NGC 6565. The large scale appearance of the Balmer ra- 
tio is the same in all the six spectra: there is a minimum 
(3.6-3.8) in the central regions, it increases outwards (up 



to 4.4-4.6 just beyond the Ha peaks), then gradually de- 
creases to 3.7-4.1 (the last statement is uncertain, because 
of the low signal-to-noise in the outer regions). 

So far, Ha/H/3 variations were reported for a few PNe 
by comparing direct, narrow-band imagery, and ascribed 
to the presence of absorbing dust within (or around) the 
objects (for example, in NGC 6302 by Bohigas 1994, in 
NGC 6445 by Cuesta & Phillips 1999, in NGC 2440 by 
Cnesta & PhiUips 2000, and in NGC 6781 by Mavro- 
matakis et al. 2001). 

The advantage of a "spectral map" is evident: it 
shows the de-projected slice of nebula covered by the slit, 
whereas the "direct images superposition" refers to the 
PN projection on the sky. 

In general, two main factors affect the Ha/H/J ratio: 

a ) the amount of dust along the line of sight. We can 
write: 

c(H/3)tot = c(H/3)i„terstellar + c(H/3)iocal (1) 

where c(H/?)iocai is null for optically thin, density 
bounded PNe. But NGC 6565 is decidedly optically 
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Fig. 8. Isophotal contours of the Ha/H/? ratio superimposed to the Ha spectral image for some representative PA of 
NGC 6565. The orientation is as in Fig. 7. The isophotal contours range between 3.6 (the outermost) and 4.6, with a 
constant step of 0.2. 



thick to the UV stellar radiation, mainly in the dense, 
equatorial regions. Moreover, Gathier et al. (1986) and 
Stasinska ct al. (1992) suggest that the infrared and 
radio excesses of the nebula could bo due to a large 
amoimt of local absorbing material. Our Ha/H/? spec- 
tral maps indicate that the neutral gas, if present, 
mainly surrounds the equatorial regions; 
b ) the Tc distribution in the ionized gas: the higher the 
electron temperature, the smaller the intrinsic Ha/H/? 
ratio. For example, we have Hq/H/3=2.74, 2.85 and 
3.02 for Te=20000 K, 10000 K and 5000 K, respectively 
(assuming the case B of Baker & Menzel 1938, and 
log Ne=3.00; Brocklehurst 1971, AUer 1984, Hummer 
& Storey 1987). 

A radial variation of Te within NGC 6565 is expected, 
due to the complexity of the ionization structure. 

Note that a large amount of dust in the equato- 
rial regions of the nebula should produce spatial and 
wavelength-asymmetric Ha/H/9 spectral maps, whereas a 
radial Te variation should not. Unfortunately, our resolu- 
tions are inadequate (by a factor of almost 2) to test such 
a possibility. In the present case, the discrimination be- 
tween a) and b) requires a detailed analysis of the nebular 
physical conditions. 

6. The physical conditions 

The radial profiles of Te and Ne are obtained from the 
classical diagnostic line ratios (ions in p'^ configuration 
for Ne and in p^ or p^ configurations for Tc; AUer 1984, 
Osterbrock 1989) and from the absolute Ha flux distri- 
bution (Sabbadin et al. 2000a, b; Ragazzoni et al. 2001). 
A compilation of the relevant references for the transition 
probabilities and collision strengths used in this paper to 
derive Te, Ne and the ionic and chemical composition is 
given by Hyung & Allcr (1996) and Liu et al. (2000). 

In all cases we make use of the zvpc, which is indepen- 
dent on the expansion velocity field, and corresponds to 
the equatorial, densest nebular regions of NGC 6565. The 
logical steps of our analysis are the following: 

- first we obtain Te[OIII] from A5007A/A4363A (the fine 

ratio is almost independent on Ne for Ne<10^ cm~^); 

- Ne[SII] is then derived from A6717A/A673lA (using 

Te[OIII] to take into account the weak dependence of 
the ratio on Te); 

- Te[NII] comes from A6584A/A5755A (adopting Ne[SII] 

for its weak dependence on Ne); 

- last, Ne(HQ!) is given by the Ha flux distribution in the 

zvcp (both Te[OIII] and Te[NII] are considered for the 
Te dependence of the Ha emissivity). 



The other diagnostics present in the spectra (like 
A471lA/A4740A of [ArlV], A5198A/A5200A of [NI], 
A5517A/A5537A of [ClIII] and the [OH] red quartet 
around A7325A) are too weak for an accurate spatial anal- 
ysis. 

This procedure, performed at all the six PA of 
NGC 6565, gives very similar results, as expected for an 
ellipsoid seen almost pole-on. Thus, only two representa- 
tive PA are illustrated here: PA=0° (along the apparent 
major axis) and PA=90° (along the apparent minor axis). 
The results are summarized in Fig. 9. 

6.1. TefOIII] 

In order to get the intrinsic line ratio I(A5007A) / 
I(A4363A), the observed line intensities must be de- 
reddened. We have just found that Ha/H/3 systematically 
changes over NGC 6565 and, at the moment, we cannot 
identify the main cause for the variation (Te or absorp- 
tion). Thus, we will consider two extreme cases: 

A) the Ha/H/3 ratio only depends on the amount of ab- 
sorption along the line of sight. In practice this corre- 
sponds to assume a constant Te= 10000 K; 

B) the absorption is constant over the nebula and the 
Ha/H/3 variation is due to Te. We discuss the case of 
c(H/3)=0.45; the choice of other, reasonable values of 
c(H/3) doesn't modify the conclusions here obtained. 

The resulting Te[OIII] profiles, shown in Fig. 9, middle 
panel, appear very similar in the two cases: the electron 
temperature is large (14000-15000 K) in the innermost re- 
gions, it rapidly decreases outwards down to about 10000 
K, and then it slightly increases (the last statement is 
weakened by the faintness of the [OIII] auroral emission). 
Previous Te[OIII] determinations in NGC 6565 are mean 
values and range from 9700 K (Acker et al. 1991) to 10300 
K (Aller et al. 1988; de Preitas Pacheco et al. 1992). 

The knowledge of the detailed Te[OIII] radial profile 
allows us to eliminate the dependence of the Ha/H/3 ratio 
on Te, thus obtaining c(H/3). Our results indicate that: 

- c(H/^)i„terstellar=0.45(±0.02); 

- in the innermost regions c(H/3)iocai —0. Once corrected 

for the interstellar component, the ratio Ha/H/3 is low, 
due to the large Te of the ionized gas. Te decreases 
outwards and, accordingly, Ha/H/3 increases; 

- the Ha/H/3 rise beyond the Ha intensity peak is 

mainly caused by the presence of neutral absorb- 
ing material around the ionized nebula. The max- 
imum of c(H/3)iocai, 0.12(±0.03), occurs at about 
1.2(±0.1)xrH+ (rH+ being the radius of the Ha peak). 
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Fig. 9. Radial physical conditions and other parameters in NGC 6565 at PA=0 ° (left) and 90° (right). Top panel: 
the observed Ha/H/3 profile in the zvpc (left ordinate scale; continuous hne) and the adopted c(H/3)tot profile (right 
ordinate scale; dashed line). Middle panel: Te[OIII] for cases A (i.e. c(H/3) variable across the nebula; continuous line) 
and B (c(H/3)=constant=0.45; long-dashed line) and Te[NII] (dotted line). Bottom panel: Ne[SII] (thick continuous 
line) and Ne from the Ha flux for some representative values of eiXr^gpixD (arcsec Kpc). The three radial zones 
discussed in the text are indicated in the lower right panel. 



The adopted, final c(H/3)tot profile for the two selected 
PA of NGC 6565 is shown in Fig. 9, top panel. It will be 
used in the next Sections to correct the observed intensi- 
ties through the relation: 

I(A)co 



log- 



fxcim 



(2) 



I(A)obs 

where f^ is the interstellar extinction coefficient given 
by Seaton (1979). 

Notice that the adopted c(H/3)tot profile is obtained 
under the assumption Te(H+)~Te[OIII]c:iTe[NII]; this is 
questionable in the outermost (knotty) regions, where the 
ionization drops and H"*" and H° co-exist. 

Lowering Te(almost neutral gas), c(H/3)iocai decreases 
and goes to zero for Te(almost neutral gas)~5000 K. In 
this case the overall Ha/H/? spectral distribution of Sect. 5 
is due to Te variations. On the one hand this modifies only 



moderately the de-reddened line fluxes, on the other hand 
it implies that H+ 3>H° (in other words: if the electron 
temperature of the almost neutral gas around NGC 6565 
is ~5000 K, there is no neutral gas!). We will deep the 
question in the next Sections. 



6.2. Ne[SII] 

The intensity ratio A6717A/A673lA of [SII] depends 
strongly on Ne and weakly on Te (Aller 1984; Osterbrock 
1989). The accurate knowledge of the c(H/3) profile is 
unimportant, because the two lines are very close. The 
resulting Ne[SII] distribution for the two selected PA of 
NGC 6565 is included in Fig. 9, bottom panel (thick con- 
tinuous line). Our analysis suggests that: 
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- the Nc[SII] profile is aligned to the Ha one. and both 

peak internally to the A6717A and A6731A maxima 
(which, in their turn, are lightly misaligned). This in- 
dicates that the S+ emission occurs in a recombining 
region; 

- the Nc[SII] distribution shows peaks up to 2800±200 

cm^^; 

- the [SII] density peak is anticorrelated to r^vpc; 

- the S+ emitting zone, being co-spatial to the N+ one 

(see Fig. 1), has a patchy and inhomogeneous struc- 
ture favoring the detection of the densest, i.e. bright- 
est, parts and indicating that the local filling factor, 
ei, in the external, low ionization layers of NGC 6565 
is ei<l. 

Previous Ne[SII] determinations in the nebula (mean 
values) date from Aller et al. (1988, Ne=1600 cm'^). 
Acker et al. (1991, Nc=1460 cm~^) and dc Frcitas Pachcco 
et al. (1992, Ne=1000 cm'^). Moreover, Meatheringham 
et al. (1988) obtain Ne=1550 cm'^ from A3726A/A3729A 
of [Oil]. 

6.3. Te[NII] 

The line intensity ratio A6584A/A5755A of [Nil] depends 
strongly on Te and weakly on Ne (Aller 1984; Osterbrock 
1989). Adopting the Ne[SII] distribution wc have derived 
the Te[NII] one shown in Fig. 9 (middle panel, dotted 
line). Te[OIII] and Te[NII] are only partially superim- 
posed: [OIII] is present all over NGC 6565, whereas [Nil] 
is peaked in the outermost, low ionization regions. Their 
values are similar; the electron temperature discrepancy, 
normally observed in PNe, can be ascribed to inaccuracy 
in the atomic parameters, small scale ionization fluctua- 
tion and, mainly, weakness of the auroral lines (Aller 1990; 
Gruenwald & Viegas 1995; Mathis et al. 1998). Previous 
Te[NII] determinations in NGC 6565 are mean values and 
span in the range 8900-9500 K (Aller et al. 1988; Acker 
et al. 1991; de Preitas Pacheco et al. 1992). 

6.4- Ne from the Ha flux 

Following Sabbadin et al. (2000a, b) and Ragazzoni et 
al. (2001), the accurate Ne radial distribution can be ob- 
tained from the absolute Ha flux in the zvpc. The de- 
reddcncd Ha profile must be further corrected for the in- 
strumental resolution, thermal motions and fine structure 
(for details, sec Sabbadin et al. 2000a). These corrections 
were superfluous for the diagnostics just analyzed, i.e. 
A5007A/A4363A, A6717A/A673lA and A6584A/A5755A, 
because they are forbidden lines and each ratio refers to 
the same ion and is emitted in the same nebular region. 

For a pixel of the Ha zvpc we have: 
AnD^FiHa) = 4njiiaN{H+)NeViei (3) 
where: 

- D is the nebular distance; 



- jna is the emission coefficient (case B of Baker & Menzel 

1938), interpolated from Brocklehurst (1971) values; 

- Vi is the "local volume" , i.e. the nebular volume iden- 

tified by a single pixel of the zvpc; Vi is given 
by lxbxsx(D/206265)3, where l=slit width (1.0"), 
b=pixel spatial scale (0.27") and s is the radial depth 
of the zvpc, s=rcspiAV/Vexp, with AV=pixel spectral 
resolution (1.79 km s~^), Vexp=28 km s~^ (see Table 
1) and rcspi is the (still unknown) angular radius of the 
cspl (i.e. the Ha nebular size in the radial direction); 

- ei, the "local filling factor", corresponds to the fraction 

of the local volume Vi which is actually filled by matter 
with density Ne. 

We can make the reasonable assumptions that Ne= 
1.15xN(H+) and ei=constant along the zvpc, thus obtain- 
ing: 

1-07 X 109 FjHa) 
Te-O-47 ^eixrcspixD^ 



Ne 



(4) 



Eq. (4), combined with Ne~Ne[SII], represents an im- 
portant link between the spatial and the dynamical prop- 
erties of the ionized gas (i.e. between the zvpc and the 
cspl). It can be used to derive the distance if independent 
information on the radial radius, rcspi, is available (for ex- 
ample, ei=l and rcspi=rzvpc in a homogeneous PN present- 
ing a large degree of symmetry). In the case of NGC 6565 
we have the opposite situation: the distance will be ob- 
tained in the next Section by means of the extinction- 
distance correlation; thus, Eq. (4) will be used to derive ei 
and the radial size of the cspl. 

Figure 9, bottom panel, shows the Ne distribution for 
three representative values of eixrcspixD (1.25, 2.5 and 
5.0 arcsec Kpc). Note the Ne radial asymmetry (the de- 
cline outward the maximum is steeper than the internal 
rise), and the similarity with the [SII] density profile, the 
alignment of the peaks and the close overlap for eiX rcspi 
xD=2.5 arcsec Kpc. 

6.5. Physical conditions: general considerations 

Our analysis shows that the equatorial regions of 
NGC 6565 consist of three, distinct radial zones, indicated 
in the right-bottom panel of Fig. 9: 

- the "fully ionized" nebula is the innermost one (roughly 

extending up to the Ha peak). Ne monotonically in- 
creases outwards, while Te first decreases then remains 
almost constant; 

- the "transition" zone located just beyond the electron 

density peak. Here Ne rapidly decreases to a few hun- 
dreds cm^'' and Te softly increases outwards; 

- the outermost nebula, the "halo", presenting a gentle 

gradient of Ne; nothing can be said on Tchaio- 

Note that: 

- the overall Ne profile can be obtained from the Ha flux 

in the zvpc, whereas Te is limited to the main nebula. 
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due to the weakness of the [OIII] and [Nil] auroral 
lines. Much deeper echellograms are necessary (and 
highly advisable) to determine Te up to the faintest 
regions; 

- the Ne of the halo shown in Fig. 9 (bottom panel) must 

be considered as an upper limit since we have adopted 

the same kincmatical properties for the main nebula 
and the halo. A better approximation is given by: 

jVehaio(Fig.9) iVehaio(Fig.9) 

iVehalo - ^ yex2Wi^ll/2 " 2 ^ ' 

- Te refers to the ionized gas, but there are growing 

evidences indicating the presence of a considerable 

amount of almost neutral, patchy gas (at unknown Te) 
in both the "transition" zone and the "halo" ; 

Up to now we have obtained the physical conditions in 
the equatorial regions of NGC 6565 using the zvpc. The 
linear relation between expansion velocity and nebular ra- 
dius, found in Sect. 4, suggests that the same analysis can 
be performed also in the polar regions. The main problem 
is represented by the line fluxes, which are considerably 
lower in the cspl than in the zvpc. 

Only a rough Ne profile can be derived for the (almost) 
polar regions, showing density peaks of 1000(±200) cm~'^ 
in [SII] and 500(±100) cm'^ from the Ha flux (adopting 
Te=10* K and ei=l; see subsect. 7.2). 

7. The nebular distance, radius and mass and the 
central star parameters 

7.1. Distance 

In Sect. 5 we have obtained c(II/3)intorstonar=0.45±0.02, 
corresponding to E(B-V)=0.31±0.01 (Acker 1978). The 
introduction of this value in the reddening-distance dia- 
grams presented by Gathicr ct al. (1986, weight 2), Maciel 
et al. (1986, weight 1) and de Oliveira-Abans & Faundez- 
Abans (1991, weight 1) gives a nebular distance of 2.0±0.5 
Kpc. Previous individual estimates, based on lower values 
of c(II/3), span the range 0.5 to 1.55 Kpc, whereas statis- 
tical distances arc spread between 0.9 (Gathier 1987) and 
7.7 Kpc (Kingsburgh & Barlow 1992). 

Using the scheme introduced in the previous Section 
and D=2.0 Kpc we have: 

RfuUy ionized 

= 0.035 X 0.029 pc 

Rtransition = 0.054 X 0.048 pC 

Rhalo — 0.11 pc (indicative value). 

7.2. Filling factor and fcspi 

Recalling, from sub-section 6.4, that the best match of 
Ne[SII] with Ne(Ha) occurs for eixrcspixD=2.5 arcsec 
Kpc, we immediately have eixrcspi=1.25 arcsec. Note that 
for ei ~ 1 we obtain ellipsoids compressed along the polar 
axis (rcspi=1.25 arcsec, i.e. r(equatorial)/r(polar)~2.6). 



Equally improbable arc the solutions for 1;5> ei; for ei=0.1 
we have rcspi=12.5 arcsec, that is a nebula very stretched 
along the polar axis; r(equatorial)/r(polar) ~0.26. 

A more stringent upper limit to ei, ei<0.35, results 
from the intuitive assumption r (polar) >r (equatorial), i.e. 
fcspi >i'zvpc- Moreover, by extrapolating to the polar re- 
gions the anticorrelation between Ne [SII] (peak) and r^vpc 
found for the equatorial zone (cfr. Sect. 6.2), we derive 
Tcspi— 5.0 arcsec (and ei~0.25). This is confirmed by the 
morphological analysis of the "NGC 6565-likc" PNe pro- 
jected more or less equatorial-on, contained in the cata- 
logues of Acker et al. (1992), Schwarz et al. (1992), Man- 
chado et al. (1996), and Gorny et al. (1999), e.g. NGC 6886 
(PNG 060.1-07.7; see Fig. 5), M 2-40 (PNG 024.1+ 03.8), 
K 3-92 (PNG 130.4-F03.1), M 2-53 (PNG 104.4-01.6) and 
M 1-7 (PNG 189.8+07.7). 

Thus, in the following we will adopt rcspi=5.0±0.5 
arcsec (corresponding to 0.048±0.005 pc), and ei=0.25, 
constant across the S"*" emitting region (consistent with 
the patchy structure of the low ionization layers). The 
position- velocity proportionality constant, A=Vexp/rcspi, 
results to be 5.6±0.6 km s"""^ arcsec"-'^; this agrees with 
the impression that the NTT-I-EMMI emission line struc- 
ture of NGC 6565 nicely reproduces (to within 20%) the 
tomographic map of the nebular slice covered by the slit; 
in fact, in Fig. 3, 4, 7 and 8 one arcsec corresponds to 3.11 
pixels along the dispersion (cspl) and to 3.70 pixels along 
the slit (zvpc). 

The kinematical age of the nebula, tkinOcA~^~1750 
yr, represents a lower limit to the actual age, fN6565) 
since the dynamical history of the gas is unknown (AUer 
1984; Schmidt- Voigt & Koppen 1987, and Dopita et al. 
1996). An estimate of tm565 can be obtained by assum- 
ing a nebular ejection at Vexp(superwind)<Vcxp(N6565), 
followed by a constant acceleration up to Vexp(N6565). 
Vexp(superwind) can be measured in the OH/IR sources 
(i.e. stars surrounded by thick opaque envelopes of gas 
and dust producing maser emissions), which are believed 
to be the PNe progenitors (see Habing 1996 and references 
therein). For Vcxp(superwind)=15(±5) km (David et 
al. 1993; Chengalur et al. 1993, and Sjouwerman et al. 
1998) we derive tN6565 =23 00(±300)yr. 

7.3. Nebular mass 

Given the distance, the ionized nebular mass is obtainable 
in different ways (from the H/3 flux, the radio flux and the 
observed Ne distribution; AUer 1984; Pottasch 1984; Os- 
terbrock 1989). Adopting D= 2.0 Kpc, we get Mio„=0.03 
(±0.01) Mq, that is considerably lower than the mean 
nebular mass of optically thin PNe, MpNe— 0.10-0.15 Mq 
(Pottasch 1983; Boffi & Stanghellini 1994; Dopita et al. 
1996). This confirms that NGC 6565 is ionization bounded 
and suggests the presence of a considerable amount of sur- 
rounding, neutral gas. From simple geometrical consider- 
ations we can put MpNe/Mion— (Rtot/Rion)'', where Rtot 
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is the radius of the (ionized + neutral) density peak. We 
have Rtot— l-4xRion, i-e. similar to the radius of the peak 
of c(H/3) given by the Ha/H/3 spectral maps. 

An estimate of the neutral gas density, N(H°), at Rtot 
can be obtained from E(B-V)iocai= c(H/3)iocai/l-48 (Acker 
1978), taking a "normal" gas to dust ratio (Spitzer 1978) 
and assuming a radial thickness of the absorbing layer 
comparable to Rtot - The resulting value, N(H°) ~ 3(±1) x 
10^ cm~^, agrees with all the previous indications that a 
dense, neutral envelope surrounds the ionized nebula. 

A further, independent confirmation comes from the 
detailed analysis of both the consistency and distribution 
of the dark globules superimposed on the [OIII] image of 
NGC 6565 (Fig. 1): 

- a number of intensity scans secured in the neighbour- 

hood of the globules indicates a depletion of the [OIII] 
flux up to 30%, corresponding to c(H/3)giobuics up to 
0.11 (to be compared with c(H/3)iocai=0.12±0.03 ob- 
tained in Sect. 5 from the Ha/H/? spectral maps); 

- most globules appear in the northern part of the neb- 

ula, as expected of an ellipsoid denser in the equato- 
rial regions and seen almost pole-on. Figs. 3, 4 and 7 
show that the southern polar cup of NGC 6565 is ap- 
proaching, thus the shadowing by the neutral, dense 
gas located in the equatorial regions mainly affects the 
northern portion of the image. 

All this, added to the different general morphology of 
the nebula in [OIII] and [Nil] (quite homogeneous and 
diffuse in [OIII] , inhomogeneous and patchy in [Nil] ; Fig. 
1) indicates the presence of noticeable density fluctuations 
at small spatial scale (below the seeing, ~0.7 arcsec) in the 
outermost, recombining nebular regions, which are rich of 
knots, globules and condensations. 

7.^. Central star parameters 

The accurate value of my* has been obtained from 
WFPC2 frames of NGC 6565 taken through the broad- 
band fllter F555W (Ac=5407A, bandwidth=1236 A). The 
resulting magnitude is mv*=18.88(±0.05), where the un- 
known star color is the main source of inaccuracy. Previous 
my* estimates reported in the literature span the range 
15.9 (Shaw & Kaler 1989) to 20.3 (Gathier & Pottasch 
1988). 

The HI and Hell Zanstra temperatures, TzHI and 
TzHell respectively, are derived in the classical way (AUer 
1984; Pottasch 1984; Osterbrock 1989) by comparing the 
de-reddened H/? and Hell A4686A nebular fluxes with 
(my,)o. The integrated H/3 and A4686A fluxes were ex- 
trapolated from the overall line profile at each PA, as- 
suming a circular symmetry of the nebular image. We 
obtain log F(H/?)obs = -11.20(±0.03) mWxm^^ and 
I(A4686A/H/3)obs=0.12(±0.02). The corresponding val- 
ues reported in the literature are in the range -11.27 
(photoelectric photometry; O'Dell 1962) to -11.1 (slit 



spectroscopy; Acker at al. 1991) for log F(H/?)obs! and 
0.12 (photoelectric narrow-band photometry; Kohoutek & 
Martin 1981) to 0.22 (slit spectroscopy; Acker et al. 1991) 
for I(A4686A/H/3)obs. 

We derive log(TzHl)=5.20(±0.05) and log(TzHell)=: 
5.08(±0.05), thus confirming the peculiarity already re- 
ported by Pottasch (1981), Martin (1981), Reay et al. 
(1984), Gathier & Pottasch (1988, 1989) and Jacoby & 
Kaler (1989): the Zanstra discrepancy is reversed in the 
central star of NGC 6565, being TzHI>TzHeII. 

The stellar luminosities are logL/L0(TzHI)=2.3 
(±0.3) and log L/LQ(TzHeII)=2.0(±0.3) (using the bolo- 
metric corrections by Schonberner 1981). 

The low luminosity and high temperature of the star, 
combined with the short age of the nebula, suggest that 
the stellar mass, M*, is larger than the average value 
(-0.60 Mq) of the PNe nuclei (Schonberner 1979, 1981, 
1983; Vassiliadis & Wood 1993, 1994; Blocker 1995). The 
quantiflcation of through the comparison with theo- 
retical evolutionary tracks reported in the literature gives 
no solutions for the He-burning post-AGB stars (these 
models do not predict a fast luminosity decline), whereas 
for the H-burning post-AGB stars different values are ob- 
tained: 

M,>0.70 Mq (Schonberner 1981, 1983, Iben 1984), 
Mh,>0.90 Mq (Wood & Faulkner 1986, Vassiliadis & 

Wood 1994), 

M*~0.65 Mq (Blocker & Schonberner 1990, Blocker 
1995), 

These discrepancies reflect the different choices by the 
authors of: 

- the semi-empirical AGB mass-loss law, controlling both 

the evolution along the upper AGB and the transition 
into the PNe region, 

- the residual, hydrogen rich envelope (and its correlation 

with the thermal-pulse cycle phase), and the fast post- 
AGB mass- loss, determining the horizontal crossing of 
the HR diagram, 

- the treatment of the gravo-thermal energy release and 

neutrino energy losses, defining the fading to the white 
dwarf regime. 

Although a deep analysis of the central star evolution 
is beyond the aim of this paper, our observational data 
suggest the following, qualitative picture: the massive pro- 
genitor of NGC 6565 underwent a long AGB phase (fa- 
voring the degeneration of the core), a strong superwind 
leaving a small, hydrogen rich envelope (causing a rapid 
horizontal crossing), and a quite fast final fading towards 
the white dwarf domain (highly degenerated core and neu- 
trino losses). 

According to the critical analysis by Blocker & 

Schonberner (1990), Tylenda & Stasinska (1994) and 
Blocker (1995), we assume M^NGC 6565)~0.65 Mq. 
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Fig. 10. Absolute emission line radial profiles (erg cm ^ s ^; logarithmic scale) in the zvpc of NGC 6565 at PA=0° 
and 90°. The orientation is as in Fig. 9. The nebula has been considered at a distance D=2.0 Kpc. 



8. The radial ionization structure 

The detailed knowledge of the plasma diagnostics in the 
(fully ionized + transition) zone of NGC 6565 allows us 
to solve the equations of statistical equilibrium, thus ob- 
taining the point to point ionic concentrations from the 
line intensities (Peimbert & Torres-Peimbert 1971, Barker 
1978, Aller 1984, Osterbrock 1989). 

The absolute emission profiles (erg cm^'"* s^^) in the 
zvpc at PA=0° and PA=90° arc shown in Fig. 10. 

All the main ionic species are well represented in the 
echellograms, with the exception of 0+, whose principal 
emissions (A3726A and A3729A) fall outside our spectral 
range. Thus, we are forced to derive the 0+/H+ abun- 
dance from the A7319.92A/Ha ratio, where A7319.92A is 
the strongest line of the red 0+ quartet, whose relative 
intensities are practically constant in the density range 
here considered (De Robertis et al. 1985). Since the ra- 
tio A7319.92A/(A3726A+A3729A) is in the range 0.015 to 
0.075, mainly depending on Ne (Keenan et al. 1999), the 



0+/H+ abundance is obtained only at, or close to, the 
line peaks. Also, the unfavorable position of A3967A, at 

the extreme blue edge of the frame, makes the quantitative 
Ne++/H+ analysis quite uncertain. 

At this point some caveats are in order; 

1) because of the luminosity drop of the fast evolving 
central star, the ionization and thermal structure of 
NGC 6565 are out of equilibrium: with time the Ha 
peak recedes (the ionized nebula zooms out) and the 
recombination processes prevail. They are faster in 
the high ionized species; this explains, at least qual- 
itatively, the reversion of the Zanstra discrepancy re- 
ported in Sect. 7.4: TzHI is larger than TzHell because 
H/3 maintains a "longer memory" than A4686A of the 
stellar luminosity. 

Actually, the post-AGB evolution slows down and 
the luminosity gradient rapidly decreases while ap- 
proaching the white dwarf domain (Schonbcrner 1979, 
1981, 1983; Wood & Faulkner 1986; Vassiliadis & 
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Fig. 11. The radial ionization structure of NGC 6565 at PA=0° (top) and PA=90° (bottom). Tlnce solutions are 
shown for 0"/H+, corresponding to Te(neutral)=5000 K, 7500 K and ~ Tc(ionizcd). The orientation is as in Fig. 9. 



Wood 1993, 1994; Blocker 1995). Because of the high 
gas density we can consider NGC 6565 in quasi- 
equilibrium: the H"*" recombination rate, dH+/dt(x- 
Ne^, is large enough to ensure a short "relax- 
ation" of the main nebula (the recombination time, 
trcc(yr)=l/Q^BNe~1.2xlO'''/Nc, is a few dozen years; 
Q;B=effective recombination coefficient of hydrogen). 

2) No correction for seeing is applied to the observed emis- 
sion line profile; although this causes some inaccuracy 
for the sharpest radial distributions (like [01], [Nil] 
and [SII]), it doesn't modify the general results here 
obtained. 

3) No precise information is yet available for the elec- 
tron temperature of the patchy, almost neutral gas in 
the external part of the "transition" zone, where the 
recombination and cooling processes dominate. This 
mainly concerns the [01] A6300A emission, which is 
much more sensitive to the physical conditions than 
are the [Nil] and [Oil] lines (Williams 1973). Thus, 



we are forced to consider different scenarios. Three 
possible solutions are given for 0'^/H+, correspond- 
ing to Te(neutral)=7500 K, Te(neutral)=5000 K and 
Te (neutral) ~Te (ionized) . 

Fig. 11 shows the resulting radial ionization structure 

of NGC 6565 at PA=0° (top) and PA=90° (bottom). Only 
two elements (i.e. helium and oxygen) present an almost 
complete ionization structure. Note that hydrogen is well 
represented in the main nebula, but the lack of the neutral 
term appears problematic in the outermost regions, where 
H° abounds (see Sects. 6.1 and 7.3). 

This is confirmed by the (He+/H+ + He++/H+) ra- 
dial trend: it is approximately constant (0.105±0.005) up 
to the "transition" zone and later it increases outwards. 
Such a behavior is typical of an optically thick PN ionized 
by a high temperature central star (the He^ region ex- 
tending further than the H+ one; Hummer & Seaton 1964; 
Alexander & Balick 1997). Previous He/H abundances re- 
ported in the literature for NGC 6565 (mean values over 
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the whole nebula) span in the range 0.095 to 0.108 (AUer 
et al. 1988; Koppen et al. 1991; de Freitas Pacheco et al. 
1992). 

8.1. Ionization structure for Te(neutral)c:iTe(ionized) 

This corresponds to the common assumption of using a 
single temperature (Te[OIII] and/or Te[NII]) all over the 
nebula, even in the regions where recombination starts to 
be significant. 

In this case the value (07H+ + 0+/H+ + 0++/H+) 
is quite constant across the nebula (but in the innermost 
regions, where O"'"'^ and even higher ionization stages dom- 
inate). 0/H results 6.0(±1.0) x 10^^, to be compared with 
the values of 5.9x10-^, 8.7x10-^ and 10.9x10-'' reported 
by AUer et al. (1988), Koppen et al. (1991) and de Freitas 
Pacheco et al. (1992), respectively. 

The most striking feature of Fig. 11 is the low O^/H"*" 
abundance in the outer regions of NGC 6565. This is 
a direct consequence of the "large" value adopted for 
Te (neutral), favoring the emissivity of the forbidden line 
A6300A with respect to Ha. The 0'^/H+ under- abundance 
indicates a scarce efficiency of the charge-exchange reac- 
tion 0+ + HO^O° + H+ (Wilhams 1973; Aller 1984; Os- 
terbrock 1989), implying that H+:$>H'^ (in contradiction 
with the result just obtained from He). 

All this is likely due to the unrealistic assumption 
Te(neutral)~Te(ionized), for the following reasons: 

- the low value of 0°/H+ in the outer regions of NGC 6565 

disagrees with a number of observational evidences 
mentioned before, indicating the presence of a large 
amount of neutral, dusty and patchy gas around the 
ionized nebula; 

- H+3>H'^ means that the nebula is almost optically 

thin to the UV flux of the star, i.e. the total 
nebular mass is close to the ionized nebular mass 
(Mtoi^Mion^O.OSMa). Thus, NGC 6565 should be an 
under-massive nebula ejected and excited by a massive 
central star, contradicting the theoretical and observa- 
tional evidences of a direct relation between M„cb and 
(Schonberner 1981; Pottasch 1983; Vassiliadis & 
Wood 1994; Bofla & StangheUini 1994; Blocker 1995; 
Buckley & Schneider 1995; Dopita et al. 1996). 



8.2. Ionization structure for Te(neutral)<Te(ionized) 

Lowering Tc(neutral), 0°/H+ increases in the outermost 
part, suggesting that charge- exchange reactions dominate 
in this H°-rich region (see Fig. 11). 

The density profile of H+ has been derived all over the 
nebula in Sect. 6.4. The H*^ concentration in the "transi- 
tion" zone is difficult to estimate. An indication can be 
obtained as follows. 
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Fig. 12. Reconstruction of the (ionized-|-neutral) radial 
density profile in NGC 6565 at PA=90° East sector. The 
orientation is as in Fig. 9. Thin continuous line: the ob- 
served Ne~N(H"'") profile (from the Ha flux in the zvpc), 
for ei=1.0. Thick continuous lines: the N(Htot) distribution 
in the "transition" zone (from Eq. (8)) for three values of 
Te(neutral). Dashed and dotted thin lines: the N(Htot) 
profile in the "transition zone" and in the "halo" for some 
representative values of time elapsed from the start of the 
recombination (from Eq. (10)). 



In the general form, the total abundance (relative to 
hydrogen) of an element X with atomic number n is given 
by: 



(6) 



For oxygen, neglecting the higher ionization terms, this 
can be written as: 

o _ o° + o+ + 0++ 

H ~ H°+H+ ^ ' 

Assuming O /H=constant across the nebula, we obtain 
for each radial position: 

o° , 0+ , 0++ 

ht + ht + -h^ (8) 



H° + H+ 



H+ 



o 

H 



which provides the H*^ distribution in the "transition" 
zone of NGC 6565, since 0/H is derived from the "fully 
ionized" nebula. 

The impasse caused by the uncertainties on 0°/H+, 
due to Te(neutral), can be overcame by linking the "tran- 
sition" zone with the "halo" , the outermost, very low elec- 
tron density region which is in a recombination phase since 
the "recent" luminosity drop of the star. 

The question is: how recent? A lower limit of about 
100 years comes from the Curtis (1918) description of 
NGC 6565 ("No central star can be seen" with the 0.91m 
Crossley Reflector). Moreover, the value M*~0.65 Mq 
adopted for the stellar mass (see Sect. 7.4) suggests that 
the luminosity drop indicatively started 1000 years ago. 
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This must be considered as an upper limit, since the thin- 
thick transition did not occur at the beginning of the 
stellar drop, but during the decline (the exact moment 
depending on the nebular mass; for Mnob— 0.15 M0 the 
thin-thick transition started at log L*/L©~ 2.8, i.e. about 
400 years ago). 

More precise information comes from the observed 
density profile in the "transition zone" and in the "halo" . 
The Ne depletion rate for recombination is given by: 
dNe/dt = -aBNeN{H+) (9) 
Assuming Ne~N(H+) and Te=100GG K, and neglecting 
the recombination delay due to expansion, we obtain: 
N{Htot) = Neit)/[1 - 8.2 x lO^HNeit)] (10) 
which provides N(Htot)5 the density in the "transition 
zone" and in the "halo" at the start of the recombina- 
tion phase, once is known the present Nc(t), the electron 
density at time t (in years) elapsed from the beginning of 
the recombination process. 

This is summarized in Fig. 12 for PA=90° East sector, 
i.e. along the apparent minor axis of NGC 6565. It shows: 

- the whole Ne~N(H+) profile in the zvpc (from the 

Ha flux, for ei=1.0); according to the discussion in 
Sect. 6.5, we adopt Nehaio=Nchaio(Fig- 9) / 2, 

- the N(Htot) distribution in the "transition" zone for Te 

(neutral)~Te(ionized), 5000 K and 7500 K (from Eq. 
(8)), 

- the density profile, N(Htot), in the "transition zone" and 

in the "halo" for some representative values of t (from 
Eq. (10)). 

In spite of the heavy assumptions (in particular: ae is 
a function of Te and ei<1.0 in the real nebula). Fig. 12 
suggests that: 

- the recombination of the main nebula indicatively 

started about 400 years ago; 

- the nebular mass is at least 0.10 Mq, hence most of the 

gas is presently neutral (Mio„ = 0.03 ±0.01, Sect. 7.3); 

- the electron temperature in the external, almost neutral 

regions is lower than in the "ionized" nebula (5000 
K<Te(neutral)< 7500 K). 

Note that the chronological sequence defined by the in- 
tersection of the different curves in the "transition zone" 
underestimates t for t< 100-200 years, since we have ne- 
glected the contribution of the photoionization in the in- 
nermost regions of the "transition zone" . This mainly con- 
cerns the t=50yr curve, which considerably shifts to the 
right. 

As already noticed in Sect. 6.1, an unrealistic situation 
occurs for Te (neutral) ~5000 K: the complete Ha/H/? dis- 
tribution observed in the spectral maps is due to Te radial 
variations. This implies that c(H/?)iocai— and H+ ° 
(an evident nonsense stressing the urgency of a more so- 
phisticate photo- ionization model). 

The overall (fully ionized-ftransition-|-halo) structure 
of NGC 6565 closely recalls the recombination phase of 



Table 2. Total chemical abundances in NGC 6565 



He/H=0.105(±0.003) 

O/H=6.0(±1.0)xl0-'' 

N/H=1.5(±0.3)xl0-* 

Ne/H=1.6(±0.4)xl0"'' 

S/H=6.3(±1.0)xl0"'^ 

Ar/H=1.8(±0.4)xl0-® 



Tylenda's (1983, 1986) models: only the innermost nebula 

is ionized by the fading UV luminosity of the star, whereas 
the outermost part, unattainable by the direct radiation, 
simply " remembers" the UV flux received just before the 
luminosity drop. In this scenario, the three-zones division 
scheme proposed in Sect. 6.5 marks two distinct evolution- 
ary times of the star: the "fully ionized" nebula refers to 
the present low luminosity phase, the "halo" to its past 
glory. 

The "transition" zone is a special case; on the whole 

it is recombining, but the process has been gradual (the 
outermost regions are the first to recombine), and decel- 
erated (the luminosity gradient of the star decreases in 
time). Moreover, following Tylcnda (1986) and Marten 
& Szczerba (1997) the high ionization species recombine 
faster than the low ones, and for a given ion, the shorter 
the wavelength of the transition, the faster is the decline 
of the respective emission. 

The electron temperature behavior in the "transi- 
tion" zone reflects the rate of change in kinetic energy of 
(ions-l-atoms-l-electrons) and is the net result of two con- 
trasting processes: the recombinations favor an increase of 
Te since they tend to remove slow electrons, whereas the 
losses by forbidden lines lower Te (AUer, 1984). All this, 
added to the inhomogeneous structure of the outermost 
layers, makes the situation quite complex. Following the 
time dependent ionization calculations by Harrington & 
Marionni (1976), Tylenda (1986) and Marten & Szczerba 
(1997), the electron temperature gradually decreases in 
time. In short, the "transition" zone constitutes a spatial 
and temporal link between the outer edge the "fully ion- 
ized" nebula and the recombining "halo" . 

The contracting phase of NGC 6565 will last only a 
few hundred years: with time, the dilution factor due to 
the expansion will overcome the slower and slower stel- 
lar decline, thus leading to a gradual re-growing of the 
ionization front. 

Actually, we cannot completely rule out the possibil- 
ity that the re-growing already takes place in NGC 6565. 
According to the classical Stromgren model of nebu- 
lar ionization, the contraction-expansion equilibrium oc- 
curs for — l/3d(/nL*)/dt= 2/tdyn, i-e. it depends on the 
present rate of luminosity decline of the star. The evolu- 
tionary tracks for hydrogen burning post-AGB stars by 
Schonberner (1981), Wood & Faulkner (1986), Vassihadis 
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Table 3. Input parameters of the model nebula 



Radial density profile 


N(Htot) in top panel of Fig. 13 




(cfr. Sect. 8.2) 


Filling factor 


1 


Chemical abundances: 




C, Si, CI 


AUer et al. (1988) 


He, N, 0, Ne, S, Ar 


this paper 


other elements 


PN (CLOUDY default) 


Dust 


PN (CLOUDY dcfauh) 


Star 


blackbody with Th,= 120000 K and 




logL,/L0= 2.0 



& Wood (1994) and Blocker (1995) give contrasting re- 
sults. 

Although a clear answer needs first and second epoch 
HST imagery in the main nebular emissions, our feeling is 
that NGC 6565 is still contracting, because of the observed 
reversion of the Zanstra discrepancy (or, at least, it seems 
to be still contracting, due to the time delay in the nebular 
adjustement to the changing UV luminosity of the star). 

We conclude this Section with the total chemical abun- 
dances (Table 2). They are derived in the usual way, 
i.e. multiplying the observed ionic abundances for the 
corresponding ICFs, the correcting factors for the un- 
observed ionic stages (obtained both empirically. Barker 

1983, 1986, and from interpolation of theoretical nebular 
models. Shields et al. 1981, Aller & Czyzak 1983, Aller 

1984, Osterbrock 1989). Following the critical analysis by 
Alexander & Balick (1997), we have considered the to- 
tal line intensities (i.e. integrated over the whole spatial 
profile and the expansion velocity field). 

In general, the chemical abundances of Table 2 are in 
satisfactory agreement with the values reported in the lit- 
erature (Aller et al. 1988; Koppen et al. 1991; de Freitas 
Pacheco et al. 1992), but for the N overabundance sug- 
gested by Aller et al. (1988) which is not confirmed here. 

9. The photo-ionization model 

In order to test the results obtained for the physical condi- 
tions and the ionization structure, and to disentangle the 
controversial points of the foregoing Sections, we have ap- 
plied the photo-ionization code CLOUDY (Fcrland et al. 
1998) to a mass of gas having the same density distribu- 
tion and chemical composition of our PN, and a powering 
source similar to the central star of NGC 6565. 

The usual caveat: CLOUDY is a "steady-state" model. 
It works in the "fully ionized" part of NGC 6565, which is 
in quasi-equilibrium (sec Sect. 8), but it cannot depict the 
external, recombining regions originated by the luminosity 
drop of the central star. 

We consider only a single sector of a single PA 
(East sector of PA=90°, i.e. along the apparent minor 



axis) as representative of the whole equatorial regions of 
NGC 6565 identified by the zvpc. 

The input parameters of the model nebula are given in 
Table 3. The adopted N(Htot) profile corresponds to the 
"fully ionized" part of N(H+) at PA=90° East sector, for 
ei=l (see Fig. 9). It is arbitrarily extended to the "tran- 
sition" and the "halo" in order to obtain a nebular mass 
of about 0.30 Mq. The choice of different N(H) profiles in 
the outermost regions does not modify the resulting ion- 
ization structure, since these layers cannot be reached by 
the UV stellar flux. Concerning the ionizing source, we use 
a blackbody distribution having the Zanstra Hell temper- 
ature and luminosity of the central star of NGC 6565 (see 
Sect. 7.4). Once it is convolved for a seeing of 0.65 arcsec, 
the line emissivity (per unit filling factor) of the model 
nebula at 2.0 Kpc can be compared with the emission per 
unit volume presented in Fig. 10. The results are shown 
in Fig. 13. 

The steady-state model nebula convolved for seeing 
(third panel in Fig. 13) closely reproduces the "fully ion- 
ized" part of NGC 6565 (bottom panel) up to R~1.0x 10^'^ 
cm. An even better fit might be obtained through a mod- 
est adjustment of the chemical abundances. The discrep- 
ancies begin in the "transition" zone and increase outward 
("halo"). 

The comparison of the model vs. observed nebula of- 
fers an encouraging support to the results obtained in the 
previous Sections. In particular, Fig. 13 confirms the pe- 
culiar properties of [ArlV] at A4740A, as pointed out in 
Sect. 4: the structure of Ar+^ (I.P.=40.7 eV) mimics that 
of higher ionization species since it is present in the in- 
terval (40.7-59.8) eV, the upper limit being the ionization 
potential of Ar+^, i.e. it overlaps the He++ zone. 

Neon shows the opposite behaviour of Argon: Ne++ 
is expected to within the (40.9-63.4) eV range, whereas 
in Fig. 13 the [Nelll] emission at A3967A extends to the 
medium-low ionization region (see also Sect. 4 and Figs. 5 
and 10). We have verified on a number of model nebulae; 
the abnormal extension of Ne++ only occurs in the pres- 
ence of a hot central star. According to Williams (1973), 
Butler et al. (1980), Shields et al. (1983) and Clegg et 
al. (1986), a large number of stellar photons with ener- 
gies well above the HI ionization threshold maintains a 
high degree of ionization in the nebula right up to the 
point where W prevails, and the charge-exchange reac- 
tion Ne+3-hH° t^Ne++-HH+ populates the -"^P, ^0° and 
levels of Ne^^ . The ground state is then reached via ra- 
diative decay, and the collisionally excited lines at A3868A 
and A3967A arc emitted. 

Concerning the electron temperature, the radial be- 
havior of the model nebula (Fig. 13, top panel) closely 
reproduces the observed one (Fig. 9, middle panel) up 
to the "transition" zone. The external, fast Te drop of 
the model nebula clearly refers to a steady-state situa- 
tion, which does not take into account the presence of the 
broad recombining region in NGC 6565. 
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Fig. 13. Model nebula vs. NGC 6565 at PA=90°, East sector. Top panel: radial profiles of Te (left ordinate scale) and 
Ne, N(H+), N(H°) and of the input density N(Htot) (right ordinate scale) for the model nebula in the case ei=1.0. 
Second panel: absolute radial flux distribution (erg cm~'^) of the model nebula in the main emissions. Third panel: 
same as in the second panel, after convolution for a seeing of 0.65 arcsec, reproducing the observational conditions of 
the true nebula. Bottom panel: absolute flux distribution in the main emissions of NGC 6565 at PA=90°, East sector 
(cfr. Fig. 10). 



Notice the quantitative agreement between E(Ha) of 
the model nebula, convolved for seeing, and the observed 
one. This is a consequence of the dumpiness found in 
Sect. 7.2 (i.e. ei~0.25). In fact, the former is obtained from 



the N(Htot) profile under the assumption ei—1, whereas in 
the real nebula N(H+) = ey^xNe[SII]= 0.5xNe[SII]. 

In Fig. 13 the most striking discrepancy concerns the 

radial 0++ profile: in NGC 6565 (bottom panel) it closely 
reproduces the H+ one, as expected of an optically thick 
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nebula powered by a high temperature central star (Flower 
1969, WiUiams 1973), while in both the model nebula and 
the convolved model nebula (second and third panel, re- 
spectively) H+ extends further than O"'""'" and the outer- 
most ionized regions are ©"'""'"-depleted. 

We stress that the application of the photo-ionization 
code CLOUDY to the other PA of NGC 6565 gives results 
coinciding with the ones just presented for PA=90°, East 
sector. 

We conclude that the nebular modeling is a powerful 
instrument for interpreting the observed properties of each 
FN. In the specific case of NGC 6565, the nebula is so 
peculiar that it cannot be completely reproduced even by a 
detailed and sophisticate code like CLOUDY. An evolving 
model is required, which takes into account the nebular 
reaction to the rapidly changing UV luminosity of the star. 

NGC 6565 is not a unique case. The recombination 
phase must be quite common in PNe, almost ineluctable 
in the presence of a massive, i.e. fast evolving, central star. 
Tylenda (1986) lists a dozen well-studied objects, includ- 
ing the famous NGC 7027, NGC 7293, NGC 6853 and 
NGC 6720, and Corradi et al. (2000) suggest that a re- 
combining halo is present in NGC 2438, but many other 
candidates are contained in the Acker et al. (1992) and 
Kohoutek (2000) catalogues on the basis of their optical 
appearance, spectroscopic characteristics and star faint- 
ness. Notice that the combination of large optical exten- 
sion and high surface brightness makes the Ring nebula 
the ideal laboratory for analyzing in detail the recombi- 
nation process in PNe. 

10. The spatial structure 

10.1. Tomography 

The reconstruction of the emission line structure in the 
nebular slices covered by the slit was introduced by Sab- 
badin et al. (2000a, b) and Ragazzoni et al. (2001). In 
the case of NGC 6565 we have selected A4686A of Hell, 
A5007A of [OIII] and A6584A of [Nil] as representative of 
the high, mean and low ionization regions, respectively, 
and Ha as a marker of the whole ionized gas distribution. 

The spectral images of the forbidden lines were decon- 
volved for seeing, spectral resolution and thermal motions, 
while also fine structure was considered for the recombi- 
nation lines. They were de-projected using the expansion 
law derived in subsect. 7.2, i.e. adopting Vcxp/r=5.6 km 
s~^ arcsec ~^ . The resulting tomographic maps (quite sim- 
ilar to the spectral images shown in Figs. 4 and 7) were 
assembled with the procedure described by Ragazzoni et 
al. (2001), thus obtaining the spatial structure presented 
in the next Section. 

10.2. 3-D morpho-kinematical structure 

In order to render the 3-D structure, we adopt the method 
introduced by Ragazzoni et al. (2001): a series of images 



showing the nebula from different directions, separated by 
15°. Each couple of images constitutes a stereoscopic pair, 
illuding the reader to see NGC 6565 in 3-D. 

The opaque reconstruction in Ha is shown in Fig. 14 
for two absolute flux cuts: log E(Ha)=-19.55 erg s~^ cm~'^ 
(upper part of each panel) and -18.50 erg s~^ cm~^ (lower 
part). Since E(Ha)=h i' a 

Ne N(H+) ei, the adopted 
cuts correspond to Ne~300 cm~^ (upper) and Ne~1000 
cm-3 (lower) (for Te=10^ K and ei=l). 

The complexity of the nebular structure in the faint, 
low ionization polar cups is more evident in Fig. 15, refer- 
ring to A6584A of [Nil] at the two cuts: log E(A6584A)= 
-19.40 erg s""'^ cm~^ (upper part of each panel) and -18.40 
erg s~^ cm~'^ (lower part). In this case E(A6584A)ocNe 
N(N"'"), thus a precise Ne cut cannot be assigned: Fig. 15 
simply represents the "faintest" and the "brightest" low 
ionization regions, respectively. 

NGC 6565 in the high ionization Hell emission at 
A4686A (not shown here for reasons of space) consists of 
a compact, almost tubular structure recalling the Ha one 
(high cut) presented in Fig. 14. Moreover, the spatial form 
in [OIII] at A5007A closely reproduces the Ha one (Fig. 
14), as expected of the similarity of the spectral images 
and the radial profiles (see Sects. 2, 4 and 8). Finally, the 
distribution of the neutral gas in the transition zone of 
NGC 6565 (obtained from Eq. (8) using the approximation 
Htot/H"'"~0°/0"^"'" and assuming Te(neutral gas)=6000 K 
and Tc([OIII])=11000 K) mimics the N+ structure shown 
in Fig. 15, although peaked a bit further away. 

In short, NGC 6565 is an inhomogeneous triaxial el- 
lipsoid seen almost polc-on (major axis=10.1 arcsec pro- 
jected in (5,7), intermediate axis=7.1 arcsec in (—85,7), 
minor axis=6.0 arcsec in (5,-83)). The matter along and 
close to the major axis was swcpt-up by some accelerating 
agent (high velocity wind? ionization? magnetic fields?), 
forming two fast, asymmetric and faint polar cups pro- 
jected in (18,7; i.e. they are lightly misaligned with re- 
spect to the major axis) and extending up to about 10.0 
arcsec from the star. The optical nebula is embedded in a 
large cocoon of almost neutral gas, which is the result of 
the recombining processes occurred during the luminosity 
decline of the powering star. 

The [OIII] and [Nil] appearance of NGC 6565 seen 
from different directions is illustrated in Figs. 16. It gives 
a representative sample of the possible morphologies as- 
sumed by our nebula when changing the line of view, and 
confirms its structural similarity with NGC 6886 (sug- 
gested in Sect. 3), as well with K 3-92, M 1-7, M 1-59, 
M 1-66, M 2-40 and M 2-53 (see the imagery catalogues 
by Schwarz et al. 1992, Manchado et al. 1996, and Corny 
et al. 1999). 

Also the observed kinematics depends on the line of 
view, as shown in Table 4 (upper panel). NGC 6565 can 
be either a compact, roundish, fast expanding object or 
a bipolar, elongated, slow one with all the intermediate 
cases. The morpho-kinematical range further spreads out 
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Fig. 14. Opaque reconstruction of NGC 6565 in Ha at Ne~300 cm~^ (for ei=l, upper part of eacli panel) and Ne~ 1000 

cm~^ (bottom part), as seen from 13 directions separated by 15°. The line of view is given by (0, V'), where 9 is the 
zenith angle and "0 the azimuthal angle; the upper-right image is the rebuilt-nebula seen from the Earth, i.e. from 
(0,0). Each horizontal couple represents a "direct" stereoscopic pair, allowing the reader to have 12 3-D views of the 
nebula in as many directions, all together covering a straight angle. We repeat here the instructions given by Ragazzoni 
et al. (2001): to obtain the three-dimensional vision, look at a distant object and slowly insert the figure in the field 
of view, always maintaining your eyes parallel. Alternatively, you can use the two small dots in the upper part of the 
figure as follows: approach the page till the two dots merge (they appear out of focus); then recede very slowly, always 
maintaining the two dots superimposed, till the image appears in focus. 

Fig. 15. Same as Fig. 14, but for [Nil], showing the faintest and the brightest (upper and lower part of each panel, 
respectively) low ionization regions of NGC 6565. 

Fig. 16. [OIII] and [Nil] appearance of NGC 6565 from 6 directions separated by 30°, showing the different mor- 
phologies when changing the line of view. The upper-right panel [6, 1/^=0,0) corresponds to the rebuilt-nebula 
seen from the Earth (to be compared with Fig. 1). The scale is the same of Figs. 14 and 15. Remember that 
projection(6',V')=projection(0±18O°,'!/'±18O°). 



in the case of a "butterfly" PN, like Hb 5, Mz 3 and 
NGC 6537. 

Decreasing of the spatial resolution exacerbates the 
problem. For example, by putting NGC G565 at the dis- 
tance of the Large Magellanic Cloud (i.e. 50 Kpc; Ko- 
vacs 2000 and references therein) we obtain the results 
of Table 4, bottom panel. In this case the echellograms 
completely lose the spatial information due to the neb- 
ula compactness. 2Vexp[0III](peak separation) cannot be 
measured any longer (in Table 4 it has bcc^ii replaced by 
2Vexp[0III](FWHM)) and the morphological classifica- 
tion is quite uncertain. 

All this evidenziates the difficulties implied in the anal- 
ysis of the high dispersion spectra and questions the re- 
liability of the interpretations, evolutionary correlations 
and spatio-kinematical models based on echellograms at 
inadequate spatial and spectral resolutions. 

11. Discussion and Conclusions 

The long travel dedicated to the study of NGC 6565 is 
over. We started with the kinematical properties, passed 
through the physical conditions and arrived to the dis- 
tance and radius of the nebula and to the temperature and 
luminosity of the star. The radial ionization structure was 
later explored, using both the "classical" method and the 
photo-ionization code CLOUDY, and the spatial structure 
investigated by assembling the tomographic maps. 

In mice: NGC 6565 is a yoimg (2000 2500 yr), patchy, 
optically thick ellipsoid with extended polar cups, seen 
almost pole-on. It is in a recombination phase, because of 
the luminosity drop of the massive powering star, which 
is reaching the white dwarf domain. The stellar decline 
started about 1000 years ago, but the nebula remained 
optically thin for other 600 years before the recombination 
phase occurred (thus, at Galileo Galilei times NGC 6565 



was at high excitation, a bit larger and much brighter than 
at present). In the near future the ionization front will re- 
grow, since the dilution factor due to the expansion will 
prevail on the slower and slower luminosity decline. 

The nebula, at a distance of 2.0 Kpc, can be divided 
into three radial zones: 

- the "fully ionized" , extending up to 0.029-0.035 pc 
at the equator (~0.050 pc at the poles), 

- the "transition" , up to 0.048-0.054 pc (~0.080 pc at 
the poles), 

- the recombining "halo", detectable up to 0.110 pc. 
The ionized mass (~0.03 Mq) is only a fraction of the 

total mass (Mtot>0.15 Mq), which has been ejected by an 
equatorial enhanced superwind of 4(±2)xl0~^ yr~^ 
lasted for 4(±2) x 10^ years. 

From all the points of view NGC 6565 results to be 
an exciting laboratory which deserves further attention. 
Amongst the many and important aspects left unresolved 
by this paper we mention: 

- the neutral gas distribution and kinematics. Detailed 

radio and infrared observations are highly advisable; 
for instance, the spectroscopically resolved structure 
in the II2 emission at 2.122 /im, which is a signpost of 
the bipolar structure (Kastner et al. 1996); 

- the filling factor. We obtain ei=0.25 in the external, low 

ionization regions by comparing Ne(IIa) and Ne[SII]. 
The same analysis should be performed in the internal, 
high ionization parts using [ArlV] (A471lA/A4740A) 
and [ClIII] (A5517A/A5537A). Deeper echellograms 
are required; 

- the faint and fast (shocked?) polar cups, overlooked in 

our study mainly focused on the equatorial regions; 

- the mechanisms and physical processes forming and 

shaping a PN like NGC 6565. This point appears prob- 
lematic, mainly because of the huge number of models 
proposed for bipolar PNe. 
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Table 4. Morpho-kinematical characteristics 



NGC 6565 at 2.0 Kpc 



Projection 



2Vexp[OIII] 2Vcxp[NII] 
(peak separation) 



[km s 



[km s 



along the major axis 57(±2) 61 (±2) 

along the intormodiato axis 42(±2) 46(±2) 
along the minor axis 34(±2) 39(±2) 



2Vexp[OIII] 2Voxp[NII] 
(10% maximum intensity) 
[km s-i] [km 8"^] 



90(±4) 

50(±4) 
40(±4) 



105(±4) 

77(±4) 
62(±4) 



morphology apparent size 



elliptical 

bipolar 
bipolar 



[arcsecj 

10.0x8.5 

20.0x8.5 
20.0x10.0 



Projection 



NGC 6565 at the distance of the Large Magellanic Cloud 



2Vexp[OIII] 
(FWHM) 
[km s-^1 



2Vcxp[NII] 
(ijcak separation) 
[km s-i] 



along the major axis 57(±4) 47(±4) 

along the intermediate a^xis 44(±4) 36(±4) 
along the minor axis 36(±4) 31(±4) 



2Voxp[OIII] 2Vcxp[NII] 
(10% maximum intensity) 



[km S-' 

80(±4) 
60(±4) 
48(±4) 



[km S-' 

86(±4) 
66(±4) 
55(±4) 



morphology apparent size 

[arcsec] 

0.40x0.34 
0.80x0.34 
0.80x0.40 



They include: single AGB progenitors (+ fast wind + 
magnetic fields + photoionization), a planetary sys- 
tem, a binary system undergoing common envelope 
evolution, a close companion presenting a dense ac- 
cretion disk (for details, see Prank 1999 and Soker & 
Rappaport 2001). 

In all these models the nebular shaping is a slow and 
gradual process occurring in a time scale comparable to 
the PN life, whereas the superior spatial resolution of 
the HST imagery has recently revealed that bipolarity 
is a common feature already in proto-PNe (see Sahai 
& Traugcr 1998, and references therein). This, on the 
one hand means that some imprint agent acts in the 
late AGB and/or early post- AGB phase, driving most 
of the subsequent nebular evolution, on the other hand 
stresses the inadequacy of the current models. 

In addition to the specific analysis of our nebula, the 
aim of the present work was the search for a satisfying re- 
duction procedure allowing us to exploit the huge amount 
of information contained in the high resolution spectra. 
They essentially cover two wide fields: kinematics (includ- 
ing expansion velocity, tomography and spatial structure), 
and physical conditions (radial profile of Ne, Te, ioniza- 
tion, chemical abundances etc.). The difficulties connected 
to the determination of the different parameters are per- 
fectly synthesized in the AUer's (1994) sentence: "A neb- 
ula is a three-dimensional structure for which we obtain a 
two-dimensional projection" . With the introduction of the 
3-D methodology, the second half of the sentence becomes: 
"... for which we obtain the three-dimensional structure". 

The accuracy of the 3-D reconstruction is defined by 
the "relative" spectral and spatial resolutions (Ragazzoni 



et al. 2001). The first is given by RR=AV/Vexp, AV being 
the spectral resolution. RR mainly depends on the intru- 
mcntation, since the expansion velocity range of the PNe 
is quite sharp (indicatively: 20 km s~^<Vexp<40 km s~^). 
The relative spatial resolution, SS=d/Ad (d=angular ex- 
tent, Ad=seeing+guiding), essentially depends on the 
target, due to the large spread of apparent sizes ex- 
hibited by the PNe. In the case of a compact object 
like NGC 6565, the ESQ NTT+EMMI echellograms give 
RR~5 and SS~10; a significant increase of the spatial res- 
olution is expected for the largest nebulae of the sample 
(like NGC 7009 and NGC 6818). 

In conclusion, we believe that the crucial point for the 
physical interpretation of the PNe is the detailed knowl- 
edge of the radial density profile, now obtainable from the 
Ha flux distribution in the zero-velocity pixel column. It 
opens the possibility of creating a realistic model for each 
expanding nebula (this includes: PNe, nova and supernova 
remnants, shells around Population I Wolf-Rayet stars, 
nebulae ejected by symbiotic stars, bubbles surrounding 
early spectral type main sequence stars etc.). 

Acknowledgements. We greatly appreciated the suggestions, 
criticisms and encouragements by Romuald Tylenda, Arsen 
Hajian, Gary Ferland, Detlef Schonberner and Vincent Icke. 
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